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We use present cosmological observations and forecasts of future experiments to illustrate the 
power of large-scale structure (LSS) surveys in probing dark matter (DM) microphysics and unveiling 
potential deviations from the standard ACDM scenario. To quantify this statement, we focus on an 
extension of ACDM with DM-neutrino scattering, which leaves a distinctive imprint on the angular 
and matter power spectra. After finding that future CMB experiments (such as COrE-|-) will not 
signihcantly improve the constraints set by the Planck satellite, we show that the next generation 
of galaxy clustering surveys (such as DESI) could play a leading role in constraining alternative 
cosmologies and even have the potential to make a discovery. Typically we find that DESI would 
be an order of magnitude more sensitive to DM interactions than Planck, thus probing effects that 
until now have only been accessible via A-body simulations. 


I. INTRODUCTION 

Dark matter (DM) is required to explain the galac¬ 
tic rotation curves, lensing and virial motions of galaxy 
clusters, observed matter power spectrum and cosmic mi¬ 
crowave background (CMB) acoustic peaks. The current 
paradigm is that DM can be well-approximated by a col¬ 
lisionless fluid, consisting of weakly-interacting massive 
particles (WIMPs) and leading to a characteristic matter 
power spectrum (P(A:) oc k~^). However, direct evidence 
for WIMPs remains elusive and it is now legitimate to 
question the validity of the standard picture. 

There is a plethora of DM models in the literature that 
exhibit small deviations from ACDM fluctuations (e.g. 
Refs. [1-24]), many of which predict additional damp¬ 
ing and/or oscillations in the P(fc). While CMB experi¬ 
ments such as Planck allow one to constrain the cosmo¬ 
logical parameters with unprecedented precision [25, 26], 
extracting the P(fc) from Planck or the next-generation 
of CMB probes (such as COrE-b [27] or PIXIE [28]) will 
be limited by the large uncertainties involved in fore¬ 
ground modelling, which hinder any angular power spec¬ 
trum analysis at large Therefore, to unravel the na¬ 
ture of DM, a direct probe of the P(fc) is needed. Here 
we show that the next generation of large-scale structure 
(LSS) surveys could provide us with key information on 
the particle properties of DM, due to their extremely high 
precision. 

Galaxy clustering surveys [29-35] have already ob¬ 
served the imprint of Baryon Acoustic Oscillations 
(BAOs), a standard ruler to measure the Hubble ex¬ 
pansion rate, H(z), and the angular diameter distance, 
Da{z). Recently, the Baryon Oscillation Spectroscopic 


^ An additional difficulty is that the Ct are the result of the convo¬ 
lution of the P(fc) with a window (Bessel) function that accounts 
for the angular scale, thus preventing one from detecting small 
features in the P(A:). 


Survey (BOSS) collaboration [36] reported a separate ex¬ 
traction of H{z) and Da{z) to a precision of 1% [29]. 
Here we show that by exploiting all of the information 
contained in the shape of the full P{k) (rather than solely 
the BAO geometrical signature [37-39]), one can test the 
validity of the ACDM model at scales below a Mpc. 

For concreteness, we focus on one specific scenario, 
in which DM scatters elastically with the active neu¬ 
trinos (hereafter, i/CDM) [1-14]. Such a DM candi¬ 
date erases small-scale perturbations through collisional 
damping [1, 3] and suppresses neutrino free-streaming in 
the early universe. This leaves a unique signature in the 
angular and matter power spectra and provides us with 
a framework to quantify the potential of future LSS sur¬ 
veys to constrain DM microphysics. We exploit both the 
current publicly available galaxy power spectrum data 
(in particular, from the WiggleZ survey [35]) and the 
expected full-shape power spectrum measurements from 
the forthcoming Dark Energy Spectroscopic Instrument 
(DESI) [40]. 

The structure of the paper is as follows. In Sec. H, 
we describe the i^CDM scenario. In Sec. HI, we com¬ 
pute up-to-date constraints using both CMB data from 
Planck and full-shape LSS data from the WiggleZ sur¬ 
vey. In Sec. IV, we perform a forecast of the sensitivity 
of planned experiments such as COrE-b and DESI to the 
i^CDM framework (and any model that generates small 
deviations from ACDM). Finally, we draw our main con¬ 
clusions in Sec. V. 


II. DARK MATTER-NEUTRINO 
INTERACTIONS 

In the ^CDM scenario, DM remains in kinetic con¬ 
tact with the neutrino sector long after the chemical 
freeze-out (see Refs. [1-14] for previous related work). 
Small-scale DM perturbations are then erased as a result 
of ongoing elastic scattering through “collisional damp¬ 
ing” [1, 3], rather than slowly clustering under gravity. 
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At the same time, neutrinos cannot free-stream as effi¬ 
ciently as in ACDM and behave more like a relativistic 
perfect fluid. The main consequences are: (i) an enhance¬ 
ment of the CMB acoustic peaks and (ii) a reduction of 
small-scale power in the matter power spectrum [5-7] . 

DM-neutrino interactions do not affect the background 
equations but they do modify the evolution of the DM 
and neutrino fluctuations. They can be implemented by 
a modification of the Euler equations. In the conformal 
Newtonian gauge^, 


^DM = ~ — S' ^/i (0 dm — Qv) 


( 1 ) 


9^ = fc V' + k - Uvj - {9^ - 0dm) 

4 3 9 

= — 0,y - —kF^3 - —/i 


15 


10 


10 ' 




where 6, 9 and cr are the density, velocity and shear per¬ 
turbations respectively, refer to higher {i > 2) neu¬ 
trino moments, V” is the gravitational potential, H is the 
conformal Hubble parameter and S = (3/4) ^dm/Pi/- 
The key quantity in Eq. (1) is /i = a ctdm-i/ ''^'DMj which 
can be written in terms of the dimensionless quantity u 
defined as 


CTDM-v 


todm 

crxh 


MOO GevJ 


This variable u describes the ratio of the DM-neutrino 
elastic scattering cross section, ctdm-i/, to the DM mass, 
mDM; normalised to the Thomson cross section, crxh 
(see Ref. [2]). In our analyses, we will consider both s- 
wave (ctdm-i/.o = constant) and p-wave oc T^) 

cross sections. Eor p-wave cross sections, we can write 
u{a) = ua~^, where u is the present-day value and a is 
the cosmological scale factor, normalised to unity today. 
The larger the value of m, the greater the suppression in 
the linear matter power spectrum with respect to ACDM, 
for a given wavenumber fc, as shown in Fig. 1 (and can 
also be seen in Refs. [5-7]). 

DM interactions leave a further imprint in the galaxy 
power spectrum through damped acoustic oscillations, 
which, in general, show up at smaller scales than those 
illustrated in Fig. 1. They were first pointed out in the 
context of DM-photon interactions (in the weak-coupling 
regime) in Ref. [2] and were later observed in interactions 
with baryons [18], neutrinos [5, 6] and dark radiation [22- 
24] . They arise because the DM fluid acquires a non-zero 
pressure as a result of interactions with the thermal bath 
and are therefore similar to the photon-baryon fluid be¬ 
fore recombination. Although they cannot be observed 
using current data, they provide a characteristic signa¬ 
ture for future experiments. 


^ In analogy to the perturbation equations governing baryon- 
photon interactions, see e.g. Ref. [41]. 


Parameter 

Prior 

Q.hhP' 

0.005 ^ 0.1 


0.01 0.99 

100 ft 

50 ^ 100 

10® A, 

1^4 

ris 

0.5 1.5 

"Treio 

0.01 ^ 0.1 

log(u) (s-wave) 

-6-> 0 

log(M) (p-wave) 

-18 ->• -11 


TABLE I. Flat priors for the cosmological parameters con¬ 
sidered here. is the baryon density, is the DM 

density, h is the rednced Hubble parameter, As is the primor¬ 
dial spectrum amplitude, ria is the scalar spectral index, Treio 
is the optical depth and u is defined in Eq. (2). 


III. CURRENT CONSTRAINTS 

To assess how powerful the constraints from future LSS 
surveys can be, we first derive the limits set by current 
CMB and galaxy clustering surveys. These will then 
serve as a benchmark for our forecasts in Sec. IV. 

To perform this analysis, the modifications shown 
in Eq. (1) are implemented in the Boltzmann code 
CLASS [42] (see also Ref. [7]) and the posterior likeli¬ 
hoods are obtained using the Markov Chain Monte Carlo 
(MCMC) code Monte Python [43]. The prior ranges 
for these parameters are listed in Tab. I. Since u can 
vary by many orders of magnitude, we select a logarith¬ 
mic prior distribution for this parameter, in contrast to 
the linear priors used in Refs. [5-7]. 

For simplicity, we assume massless neutrinos^ and fix 
the effective number of neutrino species, Veff, to the stan¬ 
dard value of 3.046 [44]. We have verified that allowing 
iVeff to vary has an impact on the value of the Hubble 
parameter, iLg, but does not change the sensitivity to the 
u parameter. 

The current CMB constraints (using Planck 2013 -I- 
WMAP polarisation data [25]) are shown in Tab. H. The 
corresponding upper limits on the DM-neutrino scatter¬ 
ing cross section (at 95% CL) are 

- 6 X 10“^^ (mDM/GeV) cm^ , (3) 

if s-wave and 

^ 2 X 10“'‘° (mDM/GeV) cm^ , (4) 

if p-wave. These results are consistent with those quoted 
by the authors of Refs. [5, 6], with the caveat that they 
did not perform a full MCMC analysis. 


^ This is in contrast to Planck, whose analysis assumes two mass¬ 
less and one massive neutrino with rriu = 0.06 eV [25]. Such a 
small neutrino mass only affects the CMB through a slight shift 
in the angular diameter distance, which can be exactly compen¬ 
sated by a decrease in 100 h of ~ 0.6 [25] . 
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S-wave 

[u = const.) 

p-wave 

{u oc r®) 

Parameter 

Planck 2013 

COrE± 

Planck 2013 

COrE± 

Qhh^ 

0.0221 ± 0.0003 

0.02223 ± 0.00004 

0.0221 ± 0.0003 

0.02222 ± 0.00004 


0.120 ±0.003 

0.1199 ±0.0005 

0.119 ±0.003 

0.1197 ±0.0005 

100 h 

70.0 ± 1.2 

67.3 ±0.2 

68.0 ± 1.2 

67.3 ±0.2 

10® As 

2.20 ±0.06 

2.207 ±0.010 

2.19 ±0.06 

2.207 ±0.010 

ns 

0.961 ±0.008 

0.9656 ±0.0017 

0.961 ± 0.008 

0.9639 ±0.0019 

^reio 

0.090 ±0.015 

0.0792 ± 0.0002 

0.090 ±0.013 

0.0790 ± 0.0002 

logioM 

< -4.04 (95% CL) 

-4.33 (95% CL) 

< -13.6 (95% CL) 

< -14.6 (95% CL) 


TABLE II. Marginalised posteriors for s-wave (left) and p-wave (right) DM-neutrino scattering cross sections set by the Planck 
2013 data (+ WMAP polarisation) (see Sec. Ill) and the COrE+ forecast (see Sec. IV). Unless otherwise indicated, the errors 
represent the 68% CL. 


We now repeat the previous analysis adding LSS data 
on the full shape of the matter power spectrum. Con¬ 
cretely, we use the galaxy clustering information from 
the WiggleZ Dark Energy Survey [35] . The WiggleZ sam¬ 
ple consists of ~ 238,000 galaxies and covers a region of 
1 Gpc^ in redshift space. Our calculations have shown 
that comparable results can be obtained from the BOSS 
DRll measurements [29]. Following a similar analysis to 
Ref. [45], we construct the likelihood function as follows: 

- 2 log[L(d„)] = , (5) 

ij 


where the covariance matrix reads 

Cij = {Phalo{ki)Phalo{kj)) — {Phaloih)) {Phalo{kj)) , (6) 

and 


Az = 


Phaloiki') Phalo,w {ki ^ ^ a) 


(7) 


In Eq. (7), Phaio{ki) is the measured galaxy power spec¬ 
trum and Phaio,w{ki,'da) is the theoretical expectation 
for the set of model parameters da, listed in Tab. I. In 
turn, P\,a\o,w{ki,da) is a convolution of the computed 
galaxy power spectrum with the survey window func¬ 
tions, W{ki,kn), and is given by 


Q \ ' iki,kn)P\xa\oikrL/^sc\,da) 

-^halOjtR 5 j — / ^ ’ 3 ' ^ 

„ ®scl 


In this equation, Osci represents the scaling, which takes 
into account that the observed galaxy redshift has to be 
translated into a distance using a fiducial model. In this 
case, we use the same values as in Ref. [46]: fib = 0.049, 
flm = 0.297, h = 0.7, Us = I and as = 0.8. The scaling 
factor is given by Ref. [45, 47]: 


DA{zfH{z) 


(9) 


The theoretical galaxy power spectrum Phaio(fc, da) is re¬ 
lated to the matter power spectrum P™‘{k,da) through 
the relation 


( 10 ) 


where b is the bias, which is assumed to be constant. We 
analytically marginalise over b as in Ref. [48]: 

,2 '^ij Phalo,w{ki,da)C^j Chalo(fcj) \ 

J2ij Phalo,w{ki, da)CPPhalo,w{kj,da) 

In Fig. 1, we show the measured galaxy power spectrum, 
Phaio{k), from WiggleZ in the four redshift bins (0.1 < 
0 < 0.3, 0.3 < z < 0.5, 0.5 < z < 0.7 and 0.7 < z < 
0.9) exploited in our analyses [35]. We also depict the 
convolved power spectrum, Pha\o,w{ki,da), as defined in 
Eq. (8) , for the ACDM fiducial model of Ref. [46] and for 
two values of the u parameter {u = 10“"* and u = 10“®, 
both in the s-wave scenario). The characteristic damping 
of the P(fc) due to the interacting DM-neutrino fluid is 
clearly visible and allows us to tighten the constraints 
with respect to the previous CMB-only analysis. 

In Tab. Ill, we present the posteriors obtained using 
the combination of WiggleZ and CMB data. We perform 
two separate analyses, including data for which: (i) k < 
fcmax = 0.12 h Mpc“^ (purely linear regime) and (ii) 
k < kjnax = 0.2 h Mpc“^ (weakly non-linear regime). 

In terms of the DM-neutrino scattering cross sec¬ 
tion (at 95% CL) with /cmax = 0.12 h Mpc“^ (k^ax = 
0.2 h Mpc“^), we obtain 

- 4 X 10"^^ (rriDM/GeV) cm^ ; 

( < 2 X 10“^^ (mDM/GeV) cm^) , (12) 

for the s-wave cross section. As we shall see in the 
next section, these bounds are competitive with those 
resulting from our forecasts for the future CMB mission 
COrE-t. 

Meanwhile, for the p-wave cross section, we obtain 

- 4 X (todm/GoV) cm^ ; 

( < 8 X 10“^^ (rriDM/GeV) cm^) . (13) 

Therefore, including data in the weakly non-linear regime 
{k < 0.2 h Mpc“^) only strengthens the constraints by 
a factor of 2 (s-wave) and 1.25 (p-wave) with respect to 
those in the purely linear regime (fc < 0.12 h Mpc“^). 


Phalo{k,da)=b^ P”^{k,da) 
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s-wave (u 

= const.) 

p-wave 

(u cx T®) 

Parameter 

fcmax =0.12 h Mpc“^ 

fcmax = 0.2 h Mpc“^ 

fcmax = 0.12 h MpC“^ 

fcmax = 0.2 h MpC“^ 


0.0220 ± 0.0003 

0.0219 ± 0.0003 

0.0219 ± 0.0003 

0.0218 ± 0.0003 


0.122 ±0.002 

0.123 ±0.003 

0.122 ±0.002 

0.123 ± 0.002 

100 h 

70.0 ± 1.1 

66.6 ± 1.0 

66.9 ± 1.1 

66.7 ± 1.0 

10® Aa 

2.19 ±0.05 

2.19 ±0.05 

2.19 ±0.05 

2.19 ±0.05 

ns 

0.956 ± 0.007 

0.956 ±0.006 

0.956 ± 0.007 

0.955 ±0.007 

^reio 

0.086 ±0.013 

0.086 ±0.013 

0.085 ±0.013 

0.085 ±0.013 

logio(^i) 

< -4.18 (95% CL) 

-4.57 (95% CL) 

< -13.7 (95% CL) 

< -13.9 (95% CL) 


TABLE III. Marginalised posteriors for s-wave (left) and p-wave (right) DM-neutrino scattering cross sections set by the 
combination of WiggleZ full-shape galaxy power spectrum measurements and Planck 2013 (-1- WMAP polarisation) data. 
Unless otherwise indicated, the errors represent the 68% CL. 


0.1 < z < 0.3 0.3 < z < 0.5 




0.5 < z < 0.7 0.7 < z < 0.9 



FIG. 1. The data points show the galaxy power spectrum, Phaio(fc), in the four redshift bins (0.1 < z < 0.3, 0.3 < z < 0.5, 
0.5 < z < 0.7 and 0.7 < z < 0.9) measured by the WiggleZ survey [35]. We also depict the convolved power spectrum, 
P\iaio,in{ki,'da), as defined in Eq. (8), for the ACDM fiducial model of Ref. [46] (solid black) and for two values of the u 
parameter: u = 10“^ (dotted red) and u = 10“® (dashed green) in the s-wave scenario. The vertical dashed line denotes the 
separation between the linear {k < 0.12 h Mpc“^) and non-linear {k > 0.12 h Mpc“^) regimes. 
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We note that, in the s-wave scenario, the bounds are as 
much as ~ 3.5 times tighter than those using only CMB 
measurements, showing the benefits of utilising the full 
shape of the P(A). The improvement is not as significant 
for the p-wave case because the suppression appears at 
larger scales (see e.g. Ref. [23]). 


Channel 6 AT AP 

(GHz) (arcmin) (/rK-arcmin) (/iK-arcmin) 


105 

10.0 

2.68 

4.63 

135 

7.8 

2.63 

4.55 

165 

6.4 

2.67 

4.61 

195 

5.4 

2.63 

4.54 


IV. FORECASTS FOR FUTURE 
EXPERIMENTS 


TABLE IV. COrE+ 4-year sensitivity. 6 is the Full Width at 
Half Maximum (FWHM) of the beam, AT and AP are the 
temperature and polarisation sensitivities respectively [27]. 


The CMB and LSS analyses in Sec. Ill allowed us to 
obtain current constraints on the DM-neutrino elastic 
scattering cross section. We will now assess the power of 
future experiments in (i) constraining DM microphysics 
and (ii) detecting small deviations from the ACDM mat¬ 
ter power spectrum in the weakly non-linear regime. 
These two analyses require slightly different method¬ 
ologies. In the first case, we construct a mock cata¬ 
logue based on the ACDM cosmology and compute the 
strongest possible upper limit on the u parameter using 
the expected sensitivity of future experiments. In the sec¬ 
ond case, the mock data assumes small but non-negligible 
DM-neutrino interactions in order to assess our ability 
to detect them and more generally, reconstruct possible 
deviations from ACDM. In both cases, we use projected 
sensitivities. 

As in the previous section, we first consider CMB ob¬ 
servables only and then include data from LSS surveys. 
We focus on two planned experiments: (i) COrE-h [27], 
a CMB space mission currently proposed for the 2015- 
2025 ESA call, and (ii) DESI [40], a multiplexed fibre- 
fed spectrograph to detect galaxies and quasars up to 
redshift z ~ 2, that is expected to run in the 2018-2022 
timeframe. 


A. COrE+ 


We first produce full mock CMB data sets (tempera¬ 
ture and E-polarisation, plus lensing). We then compute 
the fiducial angular power spectra, Cf, using the best-fit 
cosmology reported by the Planck 2015 final mission, in¬ 
cluding the TT, TE and EE spectra [26]. To these C^, we 
add a noise component consistent with each COrE-b 
channel specification and given by 


V/*' = Sjj tT^cr'^ exp 


£{i + l) 


02 

81n2 


(14) 


where correspond to the temperature or polarisation 
errors (i.e. I,JG {T,E}). The expected temperature and 
polarisation sensitivities are given in Tab. IV. 

Following Ref. [49], the effective is given by 


Xeff('*^a) — + l)/sky 

e 


D 



(15) 


where I? is a certain function of the noised power spec¬ 
tra (see Eq. (3.4) in Ref. [49]) and \C\ and \C\ repre¬ 
sent the determinants of the theoretical and observed 
covariance matrices respectively. Finally, /sky repre¬ 
sents the observed fraction of the sky (in practice, it 
weights the correlations between multipoles when the 
map does not cover the full sky). For this analysis, we 
use /sky = 0.7 [27]. 

The third step in our analysis is to compute a Gaus¬ 
sian likelihood around our fiducial spectra, using CLASS, 
tuned to obtain a 0.01% precision on the (as in 
Ref. [49], according to Eq. (15) and with a noise given by 
Eq. (14)). Then, assuming a 4-year sensitivity and using 
Monte Python to sample the parameter space with the 
priors given in Tab. I, we can predict the sensitivity of 
COrE-l- to the ACDM parameters. Note that we only 
consider the TT, TE and EE observables. For simplicity, 
we neglect tensor modes (i.e. BT, BE and BB) as they 
have currently not been observed ]50]. 

The results are presented in Tab. II. We infer that the 
future sensitivity of COrE-b to a DM-neutrino coupling 
would be (at 95% CL) 

^ 3 X 10“^^ (moM/CeV) cm^ , (16) 

if s-wave and 

^ 2 X 10-41 (moM/CeV) cm^ , (17) 

if p-wave. 

While we find that the standard cosmological parame¬ 
ters will be measured to much higher precision than with 
Planck, there is only a modest gain in sensitivity to the 
DM-neutrino cross section. Furthermore, these limits 
are slightly worse than those obtained after combining 
Planck with current LSS data in the weakly non-linear 
regime. From these results, we expect detection with 
COrE-b to be possible for u ~ lO-^. 

To assess the power of COrE-b to detect and recon¬ 
struct the izCDM cosmology or similar deviations to 
ACDM, we also produce mock data sets with u = IQ-^ 
and u = 10“® as fiducial models (for s-wave interac¬ 
tions). We then attempt to reconstruct these models 
by means of the usual MCMC method. The u = 10“® 
case is presented in Fig. 2 (and similarly for p-wave with 
u = 10-14). With COrE-l—like CMB data, one may re¬ 
construct a universe with u = 10-4 with a 40% Icr error. 
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However, the u = 10“® case would provide us with CMB 
information entirely consistent with m = 0, in agreement 
with Eq. (16). Therefore, u> 5 x 10“® is the best sensi¬ 
tivity that one could achieve with CMB experiments in 
the near future. 


B. DESI 


The DESI survey [40] is expected to provide a wealth of 
information on the matter distribution (i.e. the P(fc)) in 
the Universe at relatively small scales and up to redshift 
2 ~ 2. To forecast the ability of DESI to discover new 
physics, we first compute the expected errors from the 
DESI instrument, following a Fisher matrix approach, 
which is the usual method used to forecast galaxy survey 
experiments^’^. 

The Fisher matrix is defined as the expectation value 
of the second derivative of the likelihood surface around 
its maximum. As long as the posterior distribution for 
the parameters is well approximated by a multivariate 
Gaussian function, its elements are given by [51-53] 


Faff 



■ 1^' 

d-d a, 


(18) 


where C = S+N is the total covariance, which consists of 
signal S and noise N terms. Once more, we take a fiducial 
cosmology defined by the parameters that best fit the 
Planck 2015 TT, TE, EE -|- lowP data [54] in the presence 
of DM-neutrino interactions with u = 10“® in the s-wave 
scenario and u = 10“^'^ in the p-wave scenario®. 

Assuming a Gaussian likelihood for the DESI band 
powers, the Fisher matrix can be written as: 


F. 


LSS 





91nPgg(fc) 91nPgg(fc) 

dda 


Ves{k) 


dk 

2(27r)3 


(19) 


LL. - 


2'Kk^dkdpL 

2(27r)3 


where Ves is the effective volume of the survey and given 
by 


(/c, /i) 


nPlk, /i) 
nP{k, n) + 1 


^survey 


( 20 ) 


where ^ is the cosine of the angle between the vector 
mode (fc) and the vector along the line of sight, and n 
is the galaxy number density (which is assumed to be 
constant throughout each of the redshift bins). 


^ http;//desi.Ibl.gov/ 

® http://sci.esa.liit/euclid/ 

® See Ref. [55] for more details on the Fisher matrix formalism for 
galaxy redshift surveys such as DESI. 


To perform the analysis, we divide the data in redshift 
bins of width Az = 0.1 and cut the small-scale data at 
k = 0.25 h Mpc“^ to avoid the highly non-linear regime. 
The lowest wavenumber (i.e. the largest scale), fcmin) is 
chosen to be greater than 2tt j^ where AU repre¬ 
sents the volume of the redshift shell. We note that using 
data in the non-linear regime would require numerical 
simulations of this model. This has been performed for 
specific cases in Refs. [10, 13]. As we will discuss later, 
constraints using this method are competitive with our 
DESI forecast. 

The real-space linear DM power spectrum, Pdm, is re¬ 
lated to the linear redshift-space galaxy power spectrum, 
-Pgg. by 


Pgg(fc) = PDM(fc) (6 + /? m")' , (21) 

where b is the bias relating galaxy to DM overdensities 
in real space (as in Eq. (10)) and /3 is the linear growth 
factor. 

DESI is expected to cover 14,000 deg^ of the sky in the 
redshift range 0.15 < z < 1.85. We use the values of the 
bias given in Ref. [56] for the three types of DESI trac¬ 
ers, namely bBLGiz)D{z) = 0.84 for the Emission Line 
Galaxies (ELGs), 6 lrg(-z).D(z) = 1.7 for the Luminous 
Red Galaxies (LRGs) and bQso{z)D{z) = 1.2 for the high 
redshift quasars (QSOs). Here, D{z) is the normalised 
growth factor and both the bias and the growth factor 
are assumed to vary in each redshift bin accordingly to 
these expressions. To combine the Fisher matrices from 
the three DESI tracers, we use the multi-tracer technique 
of Ref. [57]. 

For the s-wave scenario, we obtain a Icr error on the u 
parameter of 


^y(DESI) _ 3 7 ^ 7Q-6 ^ ('22) 

for the fiducial value of u = 10“®. For p-wave, we obtain: 

^y(DESi) _ 4 4 ^ 4Q-15 ^ (23) 

for the fiducial value of m = 10“^"^. Grucially, DESI will 
ensure a 2.5cr detection of DM-neutrino interactions if 
the strength of such a coupling is u ~ 10“® (or a ~ 2cr 
detection for u ~ 10“^^ if p-wave). 

Our results are summarised in Fig. 2, where the DESI 
68% and 95% GL allowed regions in the (Ddm^^jw) plane 
are shown (assuming the Planck 2015 fiducial cosmology 
plus an interaction strength of rt = 10“® if s-wave and 
u = 10“^^ if p-wave), along with the current constraints 
and the GOrE-l- reconstruction. One can clearly see the 
improvement in the extraction of a DM-neutrino cou¬ 
pling that will be provided by the next-generation LSS 
surveys. This analysis indicates that planned galaxy clus¬ 
tering surveys will provide an extremely powerful tool 
(competitive or even better than future GMB experi¬ 
ments) to test the fundamental properties of DM. 

Since the main impact of z/GDM is the damping of 
structure on small scales, one of the largest effects will 
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FIG. 2. The 68% and 95% CL allowed regions in the (rioMh^jM) plane for the s-wave (left) and p-wave (right) scenarios. Blue; 
current constraints from the combination of WiggleZ and Planck 2013 data, with fcmax = 0.12 h Mpc“^; Magenta: projected 
sensitivity of the upcoming COrE+ CMB experiment, assuming u = 10“® (or u = if p-wave); Orange: projected 

sensitivity of the DESI galaxy survey, again assuming u = 10“® (or u = 10“^"^ if p-wave), with fcmax = 0.25 h Mpc“^. 


be a reduction in the number of satellites around galax¬ 
ies such as the Milky Way. Until now, the only way to 
study interactions at these scales has been via Wbody 
simulations, which show that for DM-radiation couplings 
greater than u ~ 10“®, the number of satellites in the 
Milky Way would be much smaller than observed [10, 13]. 
Therefore, with the sensitivity of u ~ 3.7 x 10“® expected 
from DESI, we would have a handle on alternative sce¬ 
narios to ACDM that modify our cosmic neighbourhood, 
independently of the assumptions that go into computa¬ 
tionally expensive simulations. 


V. CONCLUSION 

Cosmology provides a promising tool to measure the 
particle properties of dark matter (DM). A DM coupling 
to visible or dark radiation (including neutrinos, axions, 
dark photons or any other light uncharged particle) can 
lead to strong departures from the standard ACDM cos¬ 
mology and produce visible signatures for CMB exper¬ 
iments and LSS surveys. In the specific case of DM- 
neutrino scattering, one expects an enhancement of the 
CMB acoustic peaks due to the fact that DM is strongly 
coupled to neutrinos and vice versa, which delays the 
neutrino free-streaming epoch and alters DM clustering 
with respect to the standard ACDM picture. However, 
the largest impact is imprinted as a damping in the mat¬ 
ter power spectrum, surveyed by large-scale structure 
(LSS) galaxy surveys. 

In this study, we have looked for the optimal method to 
measure such small departures from the ACDM scenario. 
As cosmological measurements may constitute the only 
tool available to detect such effects, it is crucial to study 
the potential sensitivity of future experiments. We have 


shown that i) with current CMB measurements, one can 
probe s-wave and p-wave DM-neutrino cross sections of 
ctdm-!/,o ^ 6 X 10“^^ (mDM/GeV) cm^ and aY)M-u,2 < 
2 X 10“^° (mDivi/GeV) cm^, respectively (at 95% GL) 
and ii) by simulating a next-generation CMB experiment 
(i.e. a COrE-|—like mission) by means of a Markov Chain 
Monte Carlo analysis, one can only weakly improve on 
the current sensitivity. 

The prospects for both constraints and detection are 
far better for future galaxy surveys, such as the DESI 
or Euclid experiments. Already, current LSS data, com¬ 
bined with Planck CMB measurements, provide compet¬ 
itive constraints to those forecasted for a future CMB 
experiment such as COrE-h. Future data from the LSS 
DESI experiment alone could improve the current sen¬ 
sitivity limit by an order of magnitude, and provide an 
accurate (percent-level) measurement of the scattering 
cross section for values above that limit. Therefore, we 
have shown that galaxy clustering surveys are an excel¬ 
lent probe to detect new physics beyond ACDM. Remark¬ 
ably, future LSS experiments will be sensitive to effects 
that until now have only been accessible via Wbody sim¬ 
ulations. 


VI. ACKNOWLEDGEMENTS 

The authors would like to thank Andreu Ariho for use¬ 
ful discussions and Will Percival for precious BOSS DRll 
measurements. The authors acknowledge the use of the 
publicly available numerical codes CLASS and Monte 
Python, and the data from the Planck Legacy Archive 
and WiggleZ survey. ME is supported by Spanish grant 
FPU13/03111 of the MECD. RJW is supported by the 
STFC grant ST/K501979/1. This work was supported 







by the European Union FP7 ITN INVISIBLES (Marie 
Curie Actions, PITN-GA-2011-289442). We would like 


to thank the Universitad Autonoma de Madrid and IFT 
for hospitality while this work was completed. 


[1] C. Boehm, P. Fayet, and R. Schaeffer, Phys.Lett. B518, 
8 (2001), arXiv:astro-ph/0012504 [astro-ph]. 

[2] C. Boehm, A. Riazuelo, S. H. Hansen, and R. Schaeffer, 
Phys.Rev. D66, 083505 (2002), arXiv:astro-ph/0112522 
[astro-ph]. 

[3] C. Boehm and R. Schaeffer, Astron.Astrophys. 438, 419 
(2005), arXiv:astro-ph/0410591 [astro-ph]. 

[4] E. Bertschinger, Phys.Rev. D74, 063509 (2006), 

arXiv:astro-ph/0607319 [astro-ph]. 

[5] G. Mangano, A. Melchiorri, P. Serra, A. Cooray, 
and M. Kamionkowski, Phys.Rev. D74, 043517 (2006), 
arXiv:astro-ph/0606190 [astro-ph] . 

[6] P. Serra, F. Zalamea, A. Cooray, G. Mangano, 
and A. Melchiorri, Phys.Rev. D81, 043507 (2010), 
arXiv:0911.4411 [astro-ph.CO]. 

[7] R. J. Wilkinson, C. Boehm, and J. Lesgourgues, JCAP 
1405, on (2014), arXiv: 1401.7597 [astro-ph.CO]. 

[8] L. G. van den Aarssen, T. Bringmann, and C. Pfrommer, 
Phys.Rev.Lett. 109, 231301 (2012), arXiv: 1205.5809 
[astro-ph. CO]. 

[9] Y. Farzan and S. Palomares-Ruiz, JCAP 1406, 014 
(2014), arXiv:1401.7019 [hep-ph]. 

[10] C. Boehm, J. Schewtschenko, R. Wilkinson, C. Baugh, 
and S. Pascoli, Mon.Not.Roy.Astron.Soc. 445, L31 
(2014), arXiv:1404.7012 [astro-ph.CO]. 

[11] J. F. Cherry, A. Friedland, and 1. M. Shoemaker, (2014), 
arXiv:1411.1071 [hep-ph]. 

[12] B. Bertoni, S. Ipek, D. McKeen, and A. E. Nelson, JHEP 
1504, 170 (2015), arXiv:1412.3113 [hep-ph]. 

[13] J. Schewtschenko, R. Wilkinson, C. Baugh, C. Boehm, 
and S. Pascoli, Mon.Not.Roy.Astron.Soc. 449, 3587 
(2015), arXiv:1412.4905 [astro-ph.CO]. 

[14] J. H. Davis and J. Silk, (2015), arXiv: 1505.01843 [hep- 
ph]. 

[15] K. Sigurdson, M. Doran, A. Kurylov, R. R. Caldwell, 
and M. Kamionkowski, Phys.Rev. D70, 083501 (2004), 
arXiv:astro-ph / 0406355 [astro-ph] . 

[16] R. J. Wilkinson, J. Lesgourgues, and C. Boehm, JCAP 
1404, 026 (2014), arXiv:1309.7588 [astro-ph.CO]. 

[17] A. Dolgov, S. Dubovsky, G. Rubtsov, and I. Tkachev, 
Phys.Rev. D88, 117701 (2013), arXiv:1310.2376 [hep- 
ph]. 

[18] X.-l. Chen, S. Hannestad, and R. J. Scherrer, Phys.Rev. 
D65, 123515 (2002), arXiv:astro-ph/0202496 [astro-ph]. 

[19] C. Dvorkin, K. Blum, and M. Kamionkowski, Phys.Rev. 
D89, 023519 (2014), arXiv:1311.2937 ]astro-ph.CO]. 

[20] C.-G. Park, J.-c. Hwang, and H. Noh, Phys.Rev. D86, 
083535 (2012), arXiv:1207.3124 [astro-ph.CO]. 

[21] F.-Y. Cyr-Racine and K. Sigurdson, Phys.Rev. D87, 
103515 (2013), arXiv:1209.5752 [astro-ph.CO]. 

[22] R. Diamanti, E. Giusarma, O. Mena, M. Archidia- 
cono, and A. Melchiorri, Phys.Rev. D87, 063509 (2013), 
arXiv:1212.6007 [astro-ph.CO]. 

[23] M. Blennow, E. Fernandez-Martinez, O. Mena, J. Re¬ 
dondo, and P. Serra, JCAP 1207, 022 (2012), 
arXiv:1203.5803 [hep-ph]. 


[24] F.-Y. Cyr-Racine, R. de Putter, A. Raccanelli, 
and K. Sigurdson, Phys.Rev. D89, 063517 (2014), 
arXiv: 1310.3278 [astro-ph.CO]. 

]25] P. Ade et al. (Planck Collaboration), Astron.Astrophys. 
571, A16 (2014), arXiv: 1303.5076 [astro-ph.CO]. 

[26] R. Adam et al. (Planck Collaboration), (2015), 

arXiv: 1502.01582 [astro-ph.CO]. 

[27] F. Bouchet et al. (COrE Collaboration), (2011), 

arXiv:1102.2181 [astro-ph.CO]. 

[28] A. Kogut, D. Fixsen, D. Chuss, J. Dotson, E. Dwek, 
et ah, JCAP 1107, 025 (2011), arXiv: 1105.2044 [astro- 
ph.CO]. 

]29] L. Anderson et al. (BOSS), Mon.Not.Roy.Astron.Soc. 
441, 24 (2014), arXiv:1312.4877 [astro-ph.CO]. 

[30] L. Anderson, E. Aubourg, S. Bailey, F. Beutler, A. S. 
Bolton, et ah, Mon.Not.Roy.Astron.Soc. 439, 83 (2014), 
arXiv: 1303.4666 [astro-ph.CO]. 

[31] F. Beutler, C. Blake, M. Colless, D. H. Jones, L. Staveley- 
Smith, et al., Mon.Not.Roy.Astron.Soc. 416, 3017 
(2011), arXiv:1106.3366 [astro-ph.CO]. 

[32] C. Blake, E. Kazin, F. Beutler, T. Davis, D. Parkin¬ 
son, et ah, Mon.Not.Roy.Astron.Soc. 418, 1707 (2011), 
arXiv: 1108.2635 [astro-ph.CO]. 

[33] N. Padmanabhan, X. Xu, D. J. Eisenstein, R. Scalzo, 
A. J. Cuesta, et ah, Mon.Not.Roy.Astron.Soc. 427, 2132 
(2012), arXiv:1202.0090 [astro-ph.CO]. 

]34] W. J. Percival et al. (SDSS), Mon.Not.Roy.Astron.Soc. 
401, 2148 (2010), arXiv:0907.1660 [astro-ph.CO]. 

[35] D. Parkinson, S. Riemer-Sorensen, C. Blake, G. B. Poole, 
T. M. Davis, et al. , Phys.Rev. D86, 103518 (2012), 
arXiv:1210.2130 [astro-ph.CO]. 

[36] K. S. Dawson et al. (BOSS), Astron.J. 145, 10 (2013), 
arXiv: 1208.0022 [astro-ph.CO]. 

[37] J. Hamann, S. Hannestad, J. Lesgourgues, C. Rampf, 

and Y. Y. Wong, JCAP 1007, 022 (2010), 

arXiv: 1003.3999 [astro-ph.CO]. 

[38] E. Giusarma, R. De Putter, and O. Mena, Phys.Rev. 
D87, 043515 (2013), arXiv:1211.2154 [astro-ph.CO]. 

[39] E. Giusarma, R. de Putter, S. Ho, and O. Mena, 
Phys.Rev. D88, 063515 (2013), arXiv:1306.5544 [astro- 
ph.CO]. 

]40] M. Levi et al. (DESI), (2013), arXiv: 1308.0847 [astro- 
ph.CO]. 

[41] C.-P. Ma and E. Bertschinger, Astrophys.J. 455, 7 
(1995), arXiv:astro-ph/9506072 [astro-ph]. 

[42] J. Lesgourgues, (2011), arXiv: 1104.2932 [astro-ph.IM]. 

[43] B. Audren, J. Lesgourgues, K. Benabed, and S. Prunet, 
JCAP 1302, 001 (2013), arXiv: 1210.7183 [astro-ph.CO]. 

[44] G. Mangano, G. Miele, S. Pastor, T. Pinto, O. Pisanti, 

et ah, Nucl.Phys. B729, 221 (2005), arXiv:hep- 

ph/0506164 [hep-ph]. 

[45] B. A. Reid, W. J. Percival, D. J. Eisenstein, L. Verde, 
D. N. Spergel, et ah, Mon.Not.Roy.Astron.Soc. 404, 60 
(2010), arXiv:0907.1659 [astro-ph.CO]. 

[46] S. Riemer-Sorensen, C. Blake, D. Parkinson, T. M. 
Davis, S. Brough, et ah, Phys.Rev. D85, 081101 (2012), 
arXiv:1112.4940 [astro-ph.CO]. 



















9 


[47] M. Tegmark et al. (SDSS), Phys.Rev. D74, 123507 
(2006), arXiv:astro-ph/0608632 [astro-phj. 

[48] A. Lewis and S. Bridle, Phys.Rev. D66, 103511 (2002), 
arXiv:astro-ph/0205436 [astro-ph]. 

[49] L. Perotto, J. Lesgourgues, S. Hannestad, H. Tu, and 
Y. Y. Wong, JCAP 0610, 013 (2006), arXiv:astro- 
ph/0606227 [astro-ph]. 

[50] P. Ade et al. (Planck), (2015), arXiv:1502.02114 [astro- 
ph.COj. 

[51] M. Tegmark, A. Taylor, and A. Heavens, Astrophys.J. 
480, 22 (1997), arXiv:astro-ph/9603021 [astro-ph]. 

[52] G. Jungman, M. Kamionkowski, A. Kosowsky, and 
D. N. Spergel, Phys.Rev. D54, 1332 (1996), arXiv:astro- 


ph/9512139 [astro-ph]. 

[53] R. A. Fisher, Journal of the Royal Statistical Society 98, 
pp. 39 (1935). 

[54] P. Ade et al. (Planck Collaboration), (2015), 
arXiv: 1502.01589 [astro-ph.CO]. 

[55] H.-J. Seo and D. J. Eisenstein, Astrophys.J. 598, 720 
(2003), arXiv:astro-ph/0307460 [astro-ph]. 

[56] A. Font-Ribera, P. McDonald, N. Mostek, B. A. Reid, H.- 
J. Seo, et al. , JCAP 1405, 023 (2014), arXiv: 1308.4164 
[astro-ph. CO]. 

[57] L. R. Abramo and K. E. Leonard, 

Mon.Not.Roy.Astron.Soc. 432, 318 (2013), 

arXiv: 1302.5444 [astro-ph.CO]. 






